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Abstract. We describe an algorithm for computing the Picard-Fuchs equa- 
tion for a family of twists of a fixed elliptic surface. We then apply this 
algorithm to obtain the equation for several examples, which are coming from 
families of Kummer surfaces over Shimura curves, as studied in our previous 
work. We use this to find correspondenced between the parameter spaces of 
our families and Shimura curves. These correspondences can sometimes be 
proved rigorously. 



1. Introduction 

Let 7r : X — >• P 1 be a family of complex algebraic varieties. As s £ P 1 varies, 
the periods of the fibers X s , i.e., integrals of holomorphically varying differential 
forms against a topologically constant family of homology classes, satisfies a certain 
differential equation, known as the Picard-Fuchs equation, whose coefficients are 
rational functions. These equations and their power series solutions are interesting 
in several respects. 

Many of the previously studied examples where of families of elliptic curves with 
some extra structure. In hope of finding new applications, we decided to study the 
related families associated with Shimura curves. 

Recall that if B is a division quaternion algebra over the field Q of rational 
numbers, which is unramified at oo in the sense that B ®q K = M2 X 2(K), and if 
A4 is a maximal order in £?, the group T of norm one elements in Ai embedds in 
SL 2 (M), via the embedding of B in M 2X 2(R), a s a discrete subgroup. The quotients 
of the complex upper half plane ri by T, and more generally by finite subgroups 
r' C r, are algebraizable as moduli spaces of abelian surfaces whose endomorphism 
algebras are certain orders in Ai (so called QM abelian surfaces), together with 
some extra structures. Under mild assumptions they carry a universal family of 
such abelian surfaces. Since K3 surfaces are easier to write down, it is natural to 
consider instead the universal family of the associated Kummer surfaces obtained 
by taking the quotient under multiplication by ±1. 

In trying to write explicit equations for QM Kummer surfaces, we were led to 
study in [3] families of quadratic twists of a fixed elliptic surface (see Section 3). 
We identified 11 examples as related to QM Kummer surfaces. In each of these 
examples we have a family of varieties over P 1 for which we know that the generic 
fiber is isogenous, in a possibly complicated way, to a Kummer surface associated 
with a QM abelian surface. Consequently, there is a correspondence between the 
base spaces for these families, and a Shimura curve. In [3] we carried out a detailed 



2010 Mathematics Subject Classification. Primary: 14D07, 14J27; secondary: 14G35, 14J28. 
Key words and phrases. Picard-Fuchs equations, Elliptic surfaces, Shimura curves. 

1 



2 



AMNON BESSER AND RON LIVNE 



analysis of this correspondence, which was highly involved and required a delicate 
study of finite discriminant forms. 

The main contribution of the present work is Algorithm 2 which computes the 
Picard-Fuchs equation for a family of twists of a fixed elliptic surface, together with 
Theorem 4.3 that claims its validity. A nice feature of our alogorithm is that it 
only requires knowledge of the Picard-Fuchs equation of the elliptic surface being 
twisted. As input to our algorithm we thus need a method for computing the 
Picard-Fuchs equation for an elliptic surface. We describe an algorithm for doing 
this borrowed from a MAPLE script of F. Beukers. We prove that the algorithm 
works since we have found no record of this in the literature. 

After describing the algorithm and the proof of the main theorem, we apply the 
algorithm to the study of the families of QM Kummer surfaces described above. In 
all of the examples we expect the resulting Picard-Fuchs equation to be of degree 3 
and to further be a symmetric square (see Section 5) of a degree 2 equation. This 
turns out indeed to be the case and we list the resulting degree 2 equations. 

On the Shimra curve side, these degree 2 equations have been studied by Elkies 
in [5]. This suggests a method of discovering and verifying the correspondeces 
between the bases of the families we study and the related Shimura curve, using 
these Picard-Fuchs equations. We describe this in all the examples. We note that 
while this method falls short of a rigorous proof of the existence of a correspondence 
between the underlying moduli problems, it is far easier than the analysis carried 
out in [3]. In fact, we needed the results of the present work to exclude some 
possibilities in [3] and furthermore, one of the cases there is left unproved even 
though we know the correspondece using Picard-Fuchs techniques. 



2. The Picard-Fuchs equation 

We briefly recall the Picard-Fuchs differential equation for a family of varieties 
over a curve. For further details see for example [10]. 

Let C be a complex analytic curve and let V/C be a local system of C- vector 
spaces of dimension n over C. The analytic rank n vector bundle V := V ®c Cc 
carries a canonical connection V defined by the condition that it vanishes on sections 
of V. We fix a meromorphic vector field d/dt on C, e.g., the one associated with 
a rational parameter t if C = P 1 . Recall that d/dt induces a covariant derivative 
operator 

V dM : V -+ V . 

Let a be a meromorphic section of V. Since V has rank n, there is going to be 
a relation 

m 

J^^d/dtTa = 

with m < n and with meromorphic functions on C. We may normalize this by 
insisting that a m = 1. 

Suppose 7 e V*(U), for some open U € C, where V* is the dual of V. The 
evaluation of 7 on a, which we suggestively write as f a, is a meromorphic function 
on U and is called a period of a. Since V vanishes on sections of V it follows easily 



PICARD-FUCHS EQUATIONS OF FAMILIES OF QM ABELIAN SURFACES 



3 



that the period J a satisfies the differential equation 




which is called the Picard-Fuchs equation associated with a. 

When V comes from geometry, a bit more can be said. Suppose that tt : X — > C 
is a smooth projective family of algebraic varieties, and that V is the family of 
cohomology groups 



For some non-negative integer I. In this case, V is canonically identifies with the 
vector bundle of de Rham cohomology groups, 



and the connection V is identified with the Gauss-Manin connection on the latter 
vector bundle. If C and ir are algebraic, it follows easily that we may take a to 
be an algebraic (mcromorphic) section of V and that then the coefficients ai in 
the Picard-Fuchs equation will be rational functions on C. We will call this a 
Picard-Fuchs equation associated with the H l of the family. 

In geometric situations we may futher take a to be a section in the sub-bundle 
Tr*fl x , c and we may take 7 to be a family of Homology classes, so that the associ- 
ated period is now indeed the integral J a. 

In applications, it will always be the case that the sub-bundle Tr^fl l x ^ c will be 
of rank 1. Thus, a is determined up to a product by a rational function. The 
Picard-Fuchs equation is in some sense unique then, since we may recover easily 
the equation associated with such a product from the equation for a (see also 
Section 6 for how to remove the remaining ambiguity). 

We can also consider Picard-Fuchs equations associated with sub-local systems 
V C R l n*C provided our chosen a resides in V (g> Oc- 

We now list the local systems considered in this work. Let H be the complex 
upper half plane. We have a family of elliptic curves ir 11 : E u — > % where for r e H 
we have 



which has a constand fiber C 2 . Sec [11, § 12] for a detailed discussion. Note 
that 7r"f2^„y^ has the section dz, where z is the standard coordinate on C, whose 
associated periods are 1 and t, hence its Picard-Fuchs equation is y" = 0. 

Let F C SL2(M) be a discrete group, acting on % by fractional linear trans- 
formations. It acts on Sh in the via the standard representation of SL 2 (M) on 
C 2 . When T C SL 2 (Z) is a congruence subgroup, the quotient Xr := T\H has 
a family Ey := r\_E" of elliptic curves above it, and both are algebraizable. The 
quotient Y\Sh is a local system on Xr, isomorphic to RV^C, with tt t the induced 
projection. 

Let B be an indefinite rational quaternion algebra and let T C B x be as in 
the introduction. Let n u : A u — >• Xp be the associated universal family of abelian 
surfaces with quaternionic multiplication. Then (see [2]) the local system R 2 7r"C 



'X/C ' 



E% =C/Z<1,t) . 



We consider the resulting local system 



Sh := MV.C , 
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splits as a direct sum of a 3-dimensional local system and a system isomorphic to 
the symmetric square of Sh, Symm 2 (Sh). 

For any family ir A : A — > X of abelian surfaces, let tt s : S ~ Kummer(A) — > X be 
the associated family of Kummer surfaces. Then, the local system R 2 7rf C splits as a 
sum of a 16-dimensional trivial system and R 2 7r^C. In particular, when A = A u we 
see that R 2 7rf C splits as a sum of a 19-dimensional trivial system and Symm 2 (Sh). 
Furthermore, we have n^n 2 ^ C Symm 2 (Sh) ® Ox r - Consequently the Picard- 
Fuchs equation satisfied by the periods of a relative 2-form lo on S is going to 
be of degree 3, and will be the symmetric square of a Picard-Fuchs equation of 
degree 2 associated with the local system Sh (see Section 5 for symmetric squares 
of equations). 

3. Elliptic surfaces and their Picard-Fuchs equations 

An elliptic surface, always considered over P 1 , is a smooth and connected com- 
pact complex algebraic surface E, together with a surjective morphism ir : E — > P 1 , 
such that the generic fiber is a curve of genus 1. We will always assume that the 
fibration is relatively minimal and has a given section, denoted 0. 

For all but a finite number of points seP 1 , the fiber E s = n^ 1 (s) is an elliptic 
curve. The singular locus £ = £(-E) of the fibration is the (finite) subset of P 1 over 
which the fibers are singular (namely 7r is not everywhere smooth). Kodaira [9] 
classified all possible types of singular fibers (see also [1, Chapter V.7]). 

The generic fiber of an elliptic surface may be given by a Weierstrass equation 
of the form y 2 = f (x), where f(x) — ax 3 + bx 2 + cx + d and a, 6, c, d are rational 
functions of the parameter ion P 1 . 

Given two distinct points a and (3 in P 1 , the quadratic twist E ai p at these 
points can be described in two ways. Algebraically, if E has Weierstrass equation 
y 2 = f(x) and a and (3 are finite points, then E a ^ has equation 

(3-1) ^y 2 = f(x) . 

Analytically, E a jj can be described as follows. Take the double cover B' — > P 1 
ramified at a and j3 and let E' be the pullback surface. Now quotient E' by the 
transformation which identifies the two fibers above each fiber of E with sign —1. 

Definition 3.1. Let E — >• P 1 be an elliptic surface as above with s E S = ~E(E). 
For A e P 1 — X let E s .\ be the twisted family at s and at A. These surfaces vary in 
a family TW S (E) over the X-lme V 1 (X) - S. 

The local system R 1 7r*C over P 1 — E has dimension 2. Its dual is the homological 
invariant (tensored with C) associated by Kodaira to the elliptic surface, and we 
denote it by F. 

A Picard-Fuchs equation for the H 1 of a general elliptic sufrace E, corresponding 
to the invariant differential u — dx/y, can be computed using Algorithm 1. It is 
taken from a MAPLE script of F. Beukers (see Section 9). We failed to find it 
documented anywhere so we give a short proof that it indeed works. 

Proposition 3.2. Algortithm 1 gives the Picard-Fuchs equation for the elliptic 
surface E. 
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Algorithm 1: Computing a Picard-Fuchs equation for an elliptic surface 
Input: An elliptic surface given by a Weierstrass equation 

y 2 = ax 3 + bx 2 + cx + d, with a, b, c, d rational functions of t 
Output: The Picard Fuchs equation y" + c\y' + c^y = satisfied by the 
periods of the invariant differential 10 — dx/y 

f <- ax 3 + bx 2 + cx + d; 
■It ^ dt ' 
Jtt ^~ St' 

q <- q A x 4 + q 3 x 3 + q 2 x 2 + q x x + q ; 

e <- =Jfl + 2f ^ Cl ZlZ + C2 /2 + _ / . fe; 

C <- COEFFICIENTS (e, x) ; 

(ci, 02,90,91,92,93, 94) <- SDLVE(C = 0); 



Proof. We express y in terms of x as y = /(x) 1//2 . Applying the covariant Gauss- 
Manin differentiation with respect to t amounts to differentiating (after eliminating 
y) with respect to t. On the invariant differential lj = f{x)~ 1 l 2 dx we find 

1 -3 

Vd/d t o; = --/ 2 / t dx 

Now we may write the general differential operator of degree 2 applied to u, 

_ / 3 5 1 3 1 3 1 \ 

( 3 - 2 ) V 2 /dt cj+ Cl V d/dt w+c 2 w = ( -/" 5 / t 2 - -.r*/u - 2 C i/" 5 /« + c 2 .T 3 j cfe ■ 

For the appropriately chosen ci and C2 this will give a trivial de Rham cohomology 
class on E/P 1 . Reduction theory (see for example [8]) tells us that it is going to be 
the differential of a rational function of the form q(x)/y n and examining the poles 
at the 2-torsion points shows that one can take n = 3 and 9 a polynomial of degree 
at most 4. This is given by 

(") -(-!*■/-»•/.+*/-•)* 

To find the Picard-Fuchs equation, we equate (3.2) to (3.3), multiply by / 5 / 2 to 
clear denominators. This gives the quantity e in the algorithm. Then we simply 
solve e = 0, identically in x, expressing the c's and 9's in terms of t. □ 

4. The Picard-Fuchs equation for a family of twists 

In this section we prove our main theorem, Theorem 4.3, which describes the dif- 
ferential equation satisfied by the periods of the H 2 of the family of twists TW S (E) , 
described in Definition 3.1, of a fixed elliptic surface E. We will in fact show that 
the periods for this H 2 can be recovered from the periods for H 1 of E and the 
differential equation can be recovered solely based on a Picard-Fuchs equation for 
H 1 of E. 
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To simplify the notation, we assume that s = 0, i.e., that the twists are at 
and a varying point. Recall from the description following (3.1) that Eq,\ can be 
obtained from F as follows: One takes a double covering tt\ : B' — > P 1 which is 
ramified exactly over and A. Let d\ : B' — > B' be the deck trasformation of the 
covering. One considers the pullback tt\E and takes the quotient Tt\E/D\ where 
D\ is the map (s, e) n- (d\(s), — e), i.e., the map that identifies the hbers at s and 
d\(s) but via the map —1. The result may have singularities at the fixed points 
and A of d\ and resolving them one obtains F 0j a- We henceforth ease notation and 
write E\ for Eo t \. 

We now write a homology class T\ € H 2 (E\, C). We will obtain T\ by modifying 
a fixed homology class F e H 2 (E,C)- In fact, we take T' in Fi(P x - E, F ), where F 
is the homological invariant (see Section 3). An element of ff^P 1 — E, F) consists 
of a formal sum XX7^ X «) where ji are paths in P 1 — E and Xi is a section of F 7i , 
in such a way that the obvious boundary map vanishes. Write the one form on 
the elliptic surface F, dx/y, as a family of differential forms uj t and consider the 
function d on the path 7$ given at a point i by G,(£) = w t . 

Having hxed F' we can write a family of 2 homology classes F> € H 2 (E\, C) as 
follows. Each path 7, can be pulled back to B' . If one of the pullbacks is Si then 
the other one is d\(5i). The section Xi pull back to both of these lifts. Since we 
identify the fiber at s with the fiber at d\(s) via the map —1 and since —1 acts as 
— 1 on the first homology it is clear that T' x :— ^2(Si,Xi) + ^2(d\(8i), — Xj) descents 
to the required homology class T\ € H 2 (E\, C). 

Let us now write explicitely a period associated with this homology class. We 
first need to write a holomorphic differential two form r/\ on E\. We have a two 
form on F, 77 = u t A dt. We can write an affine model for B' , the double cover of 
P 1 ramified at and A, as s 2 = t(t — A). The form s^ 1 • has the right behavior 
with respect to deck transformations and therefore descents to the required form 
f]x on E\. Note that the choice of as a multiplier eliminates the zeros that dt 
acquires from the ramified cover. We now compute the period J r ^ r]\. Let us write 
the function obtained by evaluating An easy computation shows this is equal to 



J 7i 

Our goal is now to compute a differential equation satisfied by u. In doing so, we 
will only use the fact that the Gi satisfy the Picard-Fuchs equation for the elliptic 
family F, which we recalled in section 3, 




Dividing by 2 we get the period 



u 




(4.1) 



y" + Cl (t)y' + c 2 (t)y = 0. 



PICARD-FUCHS EQUATIONS OF FAMILIES OF QM ABELIAN SURFACES 



7 



The computation is inspired by the computation in [4, 2.10]. 

Lemma 4.1. Suppose y satisfies (4.1). Then, for a fixed X, the function z = z\ 
(t(t — X))~ 1 / 2 y satisfies the equation 

z" + a x (t)z' + p x (t)z = 

with 

2t- A 



a x {t) = ci(t) + 



Px(t) = C 2 (t) + Ci(i) 



t(t - A) 

2t - A A 2 



'2i(i-A) 4t 2 (t-A) 2 
Proof. A straightforward computation. □ 

Suppose now that we are given two rational functions p(t) = p\{t) and q(t) = 
qx(t). We have 

(pz + qz')' = p'z + (p + q')z' + qz" 

If we force the relation 

(4.2) p + q'=a x q 

then we can write 

p'z +(p + q')z' + qz" = p'z + q(z" + a x z') = p'z - q(3 x z , 

by using the differential equation for z. The relation (4.2) gives p = aq — q' so that 
p' = a'q + q'a — q" and we finally end up with the relation 

(pz + qz')' = (a'q + q'a - q" - q(3)z . 

Now, we can do the following: We have u(X) = J z x (t)dt. Since z depends on 

A only through division by \Jt — A, we easily get by differentiating n times with 
respect to A below the integral sign, 

d n u _ 1 • 3 (2n - 1) ^ / z 



d n u _ 1 ■ 3 (2n-l) v f 

[ ' ' d\ n 2™ ^ J 



Lemma 4.2. T/iere zs a choice for q such that we may expand a'q + q'a — q" — q/3 
as a polynomial in (t — A) -1 with coefficients which are rational functions of X, 

(4.4) a'q + q'a - q "- q p = J2 c n (X)(t - X)~ n 

n 

Proof. First let q be the least common multiple of the denominators of a and f3 
as rational functions of t. Then, a'q + q' a — q^ — q (3 is a polynomial in t and 
can therefore also be written as a polynomial in t — A, with coefficients which are 
rational functions of A. Suppose that this polynomial has degree m. Then, we may 
simply take q = q (t — X)~ m . □ 



We may modify the paths ji to homotopic paths making sure that \Jt — X is 
single valued on each path. Also, the sums of the monodromies of the d around 
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the paths 7^ is because T' is closed. Thus, we have, 



by (4.3), with 



/ da x dq x d 2 q x . 

{-df q > + ^T ax -^- qxPx I zxdt 



= E / E c «( A )(*- A )~ n ^ 

^? c " (A)E X(^v 



df 
df 



dA« 

2' 



£ " (A) = 1-3 (2n-l) C " (A) ' 

We therefore proved the following. 

Theorem 4.3. Let E be an elliptic surface whose periods satisfy the differential 
equation (4.1). Then, algorithm 2 computes a differential equation with polynomial 
coefficients satisfied by a non-trivial period for H 2 of the family TWq{E). 



Algorithm 2: Computing a differential equation for periods of TWo(E) 

Input: A Picard-Fuchs equation y" + c\{t)y' + c 2 (t)y = for an elliptic 
surface E 

Output: A vector c such that a Picard-Fuchs equation for the family of 
twists TWo(E) is given by c„(A)£| 

<- c 2 (t) + ci(t) 2t 2 ( * t i A A) - 4t 2 ( t_ A) 2; 

q <- LCM (DENOMINATOR (a), DENOMINATOR 

polo <- a'qo + q' a - q' Q ' - q a (3; 
m <- DEG(polo); 
q^q (t-\)- m ; 

pol <- {a'q + q'a - q" - q(3) t ^ s+x ; 
c„(A) <- COEFFICIENT (s"™ in pol); 

Cn(A) <— 1-3 (2n-l) C nW: 



5. K3 SURFACES 

In [3] we studied a particluar class of elliptic fibrations. Out of the list of elliptic 
fibrations with 4 singular fibers compiled by Herfurtner [6] , we picked out the ones 
for which the twists give K3 surfaces, generically with Picard number 19. These 
K3 surfaces are then isogenous to Kummer surfaces associted with Abelian surfaces 
whose isogeny algebra is a rational quaternion algebra. We further picked out only 
the examples in which the quaternion algebra in question is indefinite There are 
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11 examples, which we list below (Table 1, see also [3, Table 1]). The method 
for deciding which families of twists correspond to quaternion algebras and the 
method for determining the discriminant of the associated algebra are detailed 
in [3, Proposition 2.4.1 and Lemma 2.5.1]. 

As discussed in Section 2, for each of the examples above, the resulting Picard- 
Fuchs equation should be of degree 3 and should be a symmetric square of a degree 2 
equation, which is a Picard-Fuchs equation for the Shimura local system descended 
to the base. In this section we verify that this is indeed the case, and we compute 
the degree 2 equations. 

Symmetric squares of differential equations are considered, for example in [10, 
Example 6.5.2]. Given a differential equation of degree 2, y" + ay' + by = 0, one 
looks for the equations satisfied by z = y 2 . The result is 

(5.1) z'" + az" + /3z' + jz = with a = 3a, /3 = 4b + 2a 2 + a', 7 = 4ab + 2b' . 

If we are given a differential equation of degree 3, we can check if it is a symmetric 
square of one of degree 2 and find the "square root" as described in Algorithm 3: 



Algorithm 3: Taking the square root of a degree 3 differential equation 
Input: A differential equation z'" + az" + f3z' + -yz = 
Output: A differential equation y" + ay' + by = whose symmetric square 
equals the input equation, if it exists 

a <— a/3; 

b<r- (/3-2a 2 -a')/4; 
c <- 7 - 4ab - 26'; 



Not surprisingly, in all 11 examples, the resulting differential equation turn out 
to be a symmetric square of an equation of degree 2. In table 1 below we list the 
examples together with the resulting equations of degree 2 (one can recover the 
degree 3 equation using (5.1)). 

The first column is the example number, which is the same as in Table 1 in [3]. 
The second column gives the types of singular fibers for the based elliptic surface and 
their locations and the third column gives the coefficients of the degree 2 equation. 
The final column gives the expected discriminant for the associated quaternion 

algebra. In the table 7 = ( 1+ ^~ 2 ~) _ an( j j _ Il+^zQ-. Conjugates for such 
elements are over Q. 

6. SCHWARZIAN DERIVATIVES 

Our goal in the rest of this work is to compare the differential equations that we 
obtained in the previous section to those obtained by Elkies in [5]. One essential 
problem is that the equation depends in an essential way on the choice of the section 
of the de Rham bundle. Even, as is the case for us, if the choice is between different 
sections of a line bundle, it still means that the periods could be multiplied by 
an arbitrary rational function. To compare two differential equations it is best to 
compare quantities which are invariant with respect to such scaling. 

This can be done as follows for equations of degree 2. Consider the quotient 
of two independent solutions. This is invariant with respet to scaling. It depends 
of course on the choice of the two solutions, but only up to a fractional linear 
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36A 2 (-1 + A 2 ) 2 


6 


10 


^,^,oo,0 

3 3 3 8 

4 ' 4 ' 4 ' 9 


27+87A+16A 2 . -3+16A+6A 2 
" 27A+58A 2 +8A 3 > " 4A 2 (27+58A+8A 2 ) 
_ 648+1824A+3157A 2 +476A 3 +48A 4 
A 2 (27+58A+8A 2 ) 2 


10 


11 


h,h,h,iv 

oo, 0, —1, 1 

3 3 3 8 
4' 4' 4' 9 


„ _ -1 + A+4A 2 t -13-22A+27A 2 
U ~ 2(-\+\») ' " 144(-1+A) 2 (A+A 2 ) 
27+5A+64A 2 +5A 3 +27A 4 
36A 2 (-1+A 2 ) 2 


6 



Table 1 . Twists of elliptic surfaces 



transformation. Applying the Schwarzian derivative, to be recalled next, removes 
this ambiguity. Our reference for this material is [5] (see also [7]). 

Definition 6.1. The Schwarzian derivative of a function z with respect to the 
parameter £ is the function [5, (13)] 



S c (z) 



2z'z'" - 3(z") 2 



(z') 2 

where derivatives are with respect to (. 
We recall the following relevant results 
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Proposition 6.2. (1) If z\ is obtained from z by a fractional linear transfor- 
mation, then S((zi) = S^(z). 
(2) If z is the quotient of a basis of solutions to the differential equation y" + 
ay' + by = 0, derivative taken with respect to (, then the Schwarzian deriv- 
ative of z, which is independent of the choice of solutions by the first part, 
is given by [5, (17)] 

S ( (z) = Ab-a 2 -2a' . 
This gives our required invariant. 

To describe the dependency of the parameter Q, it is better to consider the 
quadratic differential 

a c (z) = S c (z)(dC) 2 

We have the formula [5, (14)] 

W = (^) sd*) + Sr,«) 

and thus 

a n (z) = a c (z) + a v (C) 

so the quadratic differential a is not independent of the parameter, but has a simple 
transformation formula with respect to changing the variable. This will allow us, 
using the sigma invariant, to determine the change of variables that will take one 
equation to another. 

To allow us to guess the required change of variables more easily, we further 
study the residue of the sigma invariant and how it behaves with respect to change 
of variables. 

An honest quadratic differential f(()(d() 2 has a well defined residue which is the 
coefficient of (~ 2 in /. with a chance of variable £ = £(77) which has order n the 
residue is multiplied by a factor of n 2 . 

Suppose now that £ = rf 1 . Then 

= 2 n(n - l)(n - 2)<' 3 f njn - l)^ 2 

= rf 2 (2(n - l)(n - 2) - 3(n - l) 2 ) 

= rf 2 (2(n 2 - 3n + 2) - 3(n 2 - 2n + 1)) 

= n- 2 (l-n 2 ) 

It follows that if the coefficient of in S^(z) is a then the coefficient of n~ 2 in 
S ( (z) is 

n 2 a + (1 -n 2 ) = 1 - n 2 (l - a) 
This leads to the following. 

Definition 6.3. The schwarzian residue of a = /(CX^C) 2 ; denoted ress cr, is 1— 
the coefficient of (~ 2 in /(C)- 

and we have 

Proposition 6.4. if ( = ((i]) is a change of variables of degree n then 

vess o- n (z) = n 2 ress o~c,{z) 
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Note that at points where the differential is holomorphic the Schwarzian residue 
is 1 and not 0. In our examples, the local system pulls back to the Shimura 
local system Sh over the upper half plane H, where it is holomorphic. Thus, 
the Schwarzian residue is always 1/n 2 for some n. We call this n the Schwarzian 
index at the point. It is, of course, just the ellipticity index of the point. 

In the following table we list for the examples we have the sigma invariant (with 
respect to the parameter A, neglecting the d\ 2 term, and the Schwarzian indexes at 
points when it is bigger than 1 (where the original hbration had a singular fiber). 



1 


7,7,oo,0 
2,2,2,6 


3(945-652A+142A 2 -60A 3 +9A 4 ) 


4A 2 (27-14A+3A 2 ) 2 


2 


2,2,2,6 


20160+42008A+41331A 2 + 17388A 3 +3888A 4 


4A 2 (72+113A+36A 2 ) 2 


3 


-10,0,00, ± 
2,2,2,6 


3(25-210A+2179A 2 +216A 3 + 16A 4 ) 


(1-8A) 2 A 2 (10+A) 2 


4 


ir,0,oo,3 
2,2,2,6 


2025+4295A+9156A 2 -1809A 3 +729A 4 


12(-3+A) 2 A 2 (5+9A) 2 


5 


(5, 5, oo, 
2,2,2,4 


3840+2072A+43A 2 +220A 3 +48A 4 


4A 2 (32+13A+4A 2 ) 2 


6 


4,l,oo,0 
2,2,2,4 


3(20-33A+28A 2 -7A 3 +A 4 ) 


4A 2 (4-5A+A 2 ) 2 


7 


^,0,oo,± 
2,2,2,4 


15(l25-675A+4244A 2 +501A 3 +45A 4 ) 
4(1-5A) 2 A 2 (25+3A) 2 


8 


=±,0,oo,l 
2,2,2,4 


3(l+3A+13A 2 -6A 3 +9A 4 ) 


4(-l+A) 2 (A+3A 2 ) 2 


9 


l,-l,oo,0 
2,2,2,3 


32+49 A 2 +27A 4 


36A 2 (-1+A 2 ) 2 


10 


=f,=±,oo,0 
2,2,2,3 


648+1824A+3157A 2 +476A 3 +48A 4 


A 2 (27+58A+8A 2 ) 2 


11 


oo, 0, —1, 1 
2,2,2,3 


27+5A+64A 2 +5A 3 +27A 4 


36A 2 (-1+A 2 ) 2 



Table 2. Sigma invariants and Schwarzian indices 



7. The results of Elkies 

In [5] Elkies computes certain differential equations associated with Shimura 
curves. While this is not stated explicitly, these are exactly the Picard Fuchs 
equations associated with the Shimura local system descended to the Shimura curve 
because the quotient of the two solutions gives the coordinate r on the upper half 
plane, just as for the Shimura local system, as in Section 2. 

The briefly list the types of Shimura curves considered. For more information 
one may consult [5] or [3] (our notation is consistent with the latter reference). 
For each discriminant D (always the product of an even number of primes) the 
Shimura curve Vd is the quotient of the upper half plane by the group T of norm 
one elements in a maximal order in a quaternion algebra of discriminant D (see 
the introduction). For each prime p\D it carries a modular involution w p and these 
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involutions commute with each other so that we can also set w n — Yi P \ n w p f° r 
n\D. We let be the quotient of Vd by the group generated by all its modular 
involutions. Finally, for a prime p which does not divide D there is a modular curve 
Vd, p , which corresponds to an additional 'T (p)" structure, This retains all of the 
previous involutions but has an additional involution w p . 

In table 3 we give, for the relevant curves, the equation that Elkies finds, the 
associated sigma invariant and the Schwarzian indices at the relevant points. The 
equations of Elkies are in non-normalized form ay" + by' + cy = 0, so one needs to 
normalize first before computing the sigma invariant. 





t(t - 2)(t - 27)y" + iQ*!=203t±2i V + { n_ r_ )y = Q 

10368-7296t+3157i 2 -119t 3 +3t 4 
4(-27+«) 2 (-2+t) 2 t 2 


27,2,00,0 
2,2,2,3 




t{16t 2 + 13* + 8)y" + (24i 2 + 13* + %' + (|i + ±)y = 

192+440 t+43t 2 + 1036t 3 +960t 4 
4t 2 (8+13t+16t 2 ) 2 


Si, 5i, 0, oo 
2,2,2,4 


V{ 5 


(t 81)(* - l)ty" + (f - 82t + f )y> + (± - ±) y = 

35t 4 -3680t 3 +244242t 2 -2449444+177147 
° ~ 36(t-81) 2 (t-l) 2 t 2 


l,81,0,oo 
2,2,2,6 



Table 3. Elkies's list of differential equations 



Here <5i is a solution to the equation 16t 2 + 13t + 8 = 0. 

For discriminant 6 Elkies does not write down the equation explicitely, though 
he gives a recepy to discover one of 4 possible equations. As he notes, however, 
since there are only 3 elliptic points, the sigma invariant is uniquely determined by 
the indexes of ellipticity. Suppose that these are at t = 0, l,oo. The most general 
form of a is 

(a b c d \ . , , 9 

and one has the condition c + d = to avoid a pole of order 3 at oo. The residues 
are a, b and a + b + d at 0, 1 and oo respectively, from which all the coefficients are 
easily determined. In the case at hand, Elkies chooses the coordinate t so that the 
indices are 2,4,6 at 0, 1 and oo respectively. This gives 

_ / 3 15 103 103 

° ~ + 16(4- l) 2 + 1441 ~ 144(t- 1) 

8. COMPARISON WITH THE RESULTS OF ELKIES 

In this section we compare Elkies's list with the list of differential equations we 
obtained in Section 5. As explained in the introduction, each of these examples is 
a family of varieties over P 1 and there is a correspondence between these P 1 and a 
Shimura curve of some (computable) discriminant. 

The above correspondece is compatible with the Shimura local system Sh. This 
implies that the correspondence is going to map the sigma invariants on Elkies's 
list to the corresponding sigma invariant of the families. Here we demonstrate how 
one can use this to guess the correct correspondence. This is not a proof that 
the correspondence is the correct one, which then needs to be established by more 
precise means [3, Section 8]. It can nevertheless be used to exclude certain possible 
correspondenced (see Subsection 8.2 in the above reference). 
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No. 10 - Corresponds to discriminant 10. The correspondence has to carry the 
special points of the fibration at A = —27/4, —1/2, oo, with respective Schwarzian 
indices 2,2,2,3 to the special points t = 27, 2,oo,0 with the same respective in- 
dices for the equation that Elkies finds for Vm/(w 5} u> 2 ). It is trivial to guess the 
relation t = — 4A. A change of variables for the sigma invariants confirms this, It 
can be proved rigorously (see [3, Subsection 8.3]) that the A- line is isomorphic to 

V U /(W2,W 5 ) 

No. 5 - Corresponds to discriminant 14. The correspondence has to carry the 
special points of the fibration at A = <5, 5, oo,0 with respective Schwarzian indices 
2,2,2,4 to the special points t = <5i,(5i,0, oo with the same respective indices for 
the equation that Elkies finds for Via/ '(W7, 102)- Since 5 is a solution of the equation 
Ax 2 + 13a; + 32 = it is easy to guess the relation t = 2/A. A change of variables 
for the sigma invariants confirms this. It can be proved rigorously (see [3, Subsec- 
tion 8.1]) that the A-line is indeed isomorphic to the Shimura curve V14/ (w2 7 W7). 

No. 3 - Corresponds to discriminant 15. The correspondence has to carry the 
special points t = 1,81,0, 00 with respective Schwarzian indices 2,2,2,6 for the 
equation that Elkies finds for Vis/(w5, W3) to the special points of the fibration at 
A = —10,0, 00, 1/8 with the same indices. This can be done with the change of 
variables A = and a change of variables for the sigma invariants confirms this. 
It can be proved rigorously (see [3, Lemma 3.8.2]). 

No. 4 - Corresponds to discriminant 10. In this case we speculated (but could not 
prove) that the A-line was the curve ^10,3/(^2, W5, w 3 ). Here we show this is consis- 
tent with the Picard-Fuchs equations. According to Elkies, the curve V10, 3/(^2, W5) 
is a degree 4 cover of V* , rational with a coordinate x such that 

t _ (-6 + 6a;) 3 

~~ (1 + xf (17- lOx + 93: 2 ) 

(this is equation (57) in [5] but the 7 there should be corrected to 17, as for example 
in the computation between equations (59) and (60)). From the expression 

6 3 , (3x 2 + 5) 2 
97T8 Wlth T = 9(x 1)3 

which is also in (57) there it is easy to see that the map x — > t sends 1, 00, —1, 5 to 
0, 0, 00, 2 with multiplicities 3, 1, 2, 2 respectively, ±y — 5/3 to 27 with multiplicity 
2, the two roots of 9x 2 — 10a; + 5 = to 2 with multiplicity 1, and the two roots 
of 9a; 2 — 10.x + 17 = to 00 with multiplicity 1. Thus, the elliptic points for 
^10,3/(^2, W5) are going to be at x = 00 with multiplicity 3 and at the roots of 
the equations 9a; 2 - 10a; + 5 = and 9x 2 - 10a; + 17 = with multiplicity 2. The 
involution W3 is given by Elkies, just after (57), to be w^(x) = ^ _ x anc [ so a 
coordinate on the quotient is given by 

/ 5\ 2 , 25 

C-»(.- 5 ) =^-10x + - 

We see that the elliptic points will map to ( — 00, —20/9,-128/9, so these will be 
elliptic of degree 6, 2, 2, and in addition the ramification point is elliptic of degree 
2. We can map C to A with the correct orders by 

128 128 

3C+ 1 ! 5- 3(9a; 2 - 10a; + 17) ' 
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This is confirmed by the matching of the sigma invariants. 

Other examples correspond to discriminant 6. Some of them are directly inter- 
related. Consider examples number 6 and 8. The special points are Ai = 4, l,oo,0 
and A 2 = —1/3, 0, oo, 1 respectively with the same indices. There is a finite number 
of ways for carying one set to the other preserving the indices, and testing each 
one using the sigma invariants we get the correct transformation Ai = 1 — 1/A 2 . It 
turns out that (see [3, Subsection 8.7]) that making this change of variable makes 
the two base elliptic fibrations isogenous. 

Consider next examples 9 and 11. The special points in both cases are oo, 0, — 1, 1 
but with indices 2,3,2,2 in example 9 and 2,2,2,3 in example 11. Testing again 
the finite number of possible transformations with the sigma invariants gives A 2 = 
(1 + Ai)/(1 — Ai). It is proved in [3, Subsection 8.8] that this again makes the two 
base fibrations isogenous. 

No. 6 - For Vq/ (w2tW^) it turns out to be slightly better to work with the 
coordinate ( = 1/(1 — t) so that the elliptic points are at ( = 0, 1 and oo with 
indices 6, 2 and 4 respectively. To get the required ellipticity behavior for the A- 
line, with elliptic points at A = 4, l,oo,0 with indices 2,2,2,4, one may consider 
a degree 3 map having ramification type (2, 1) over oo, producing indices 2 and 4, 
ramification 3 above 0, producing an additional index 2 and ramification type (1, 2) 
above 1, producing one additional index 2 and an additional non-elliptic point. 
This can be aranged by a map of the form £ = cA _1 ((A — l) 3 for the appropriate c 
for which this ramifies above 1. So c is the value for which one of the roots of the 
derivative (c(A — l) 3 — A)' = 3c(A — l) 2 — 1 is mapped to ( = 1. We have for that 
root 

(A_1) 3 = A_1 
A 3A 

hence A = —1/2 and c = 4/27. Consider the equation for A to map to ( — 1. As an 
equation on A — 1 the sum of the 3 roots should be 0, hence the third root should 
be 3, so that the additional preimage of 1 is 4. Thus, the cover we found matches 
perfectly with the A-line. Summarizing, we have 

1-1 27A 
~ ? ~ 4(A~ T) 3 

This is confirmed by the sigma invariants. 

No. 9 - Here the elliptic points are at 1, — 1, oo and with indices 2,2,2 and 3. 
In trying to relate them with the elliptic points for ( it is very easy to guess the 
relation ( = A 2 , and this is confirmed by the cr-invariants. Thus, the A-line is a 
double cover of Vq ramified above the elliptic points of order 4 and 6. This was 
used in [3, Subsection 8.2] to prove that the A-line is Vq/(wq). 

No. 7 - The elliptic points are — 25/3, 0, oo, 1/5 and with indices 2 at the first 
3 points and 4 at the last point. Here we guess that the base for the family of 
twists is isomorphic to the quotient Vq, 5/(w 2 , W3, u^)- Elkies find a coordinate x 
on ^6, 5/(^2, W3) for which the action of w$ is given by [5, (37)] Ws(x) = (42 — 
55x)/(55 + 300a;). The two fixed points of this action are 7/30 and —3/5. Thus, a 
coordinate on the quotient is provided by 

(8.1) y=((x + 3/5)/(a;-7/30)) 2 . 

The map from A^(5) to Va/(w2,ws) is given by [5, Equation 36] by 

t = (l + 3x + 6a; 2 ) 2 (l -6x + 15a; 2 ) = 1 + 27x 4 (5 + 12a; + 20a; 2 ) . 
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The relation between t and £ is 

(8 ' 2) C = 1/(1 ~ t] = 27*4(5 + 12* + 20*2) ' 

The ramification above £ = is of order 6 at infinity. The ramification over £ = oo 
is of order 4 at x = and order 1 at each of the roots of 5 + 12.x + 20a; 2 . The 
ramification over £ = 1, or t = 0, is of order 1 at each of the roots of 1 — 6x + 15a; 2 
and of order 2 at each of the roots of 1 + 3a; + 6a; 2 . Thus the elliptic points of the 
cover are of order 4 at the roots of 5 + 12a; + 20a; 2 and of order 2 at each of the roots 
of 1 — 6a; + 15a; 2 . These two pairs of points are interchanged by W5. The elliptic 
points of order 4 arc mapped to y = —9/16 while those of order 2 are mapped to 
y = —24. In addition we get elliptic points at the ramification points of the covering 
at y = and y — 00, both of order 2. Now, if we guessed correctly, there would 
be a Mobius transformation sending the 4 elliptic points to the 4 singular points 
of the elliptic surface, sending y — —9/16 to A = 1/5 . It is easy to check that the 
unique transformation of this type is A = — 25y/(3(24 + y)). Composing with (8.1) 
we get 

-(3 + 5a;) 2 
5(1 - 6a; + 15a; 2 ) ' 

This, as usual, is confirmed by Pulling back the cr-invariants. Our guess can be 
proved rigorously [3, Subsection 8.6]. 

No. 2 - The elliptic points are at —9/4,— 8/9, 00, with indices 2,2,2,6 respec- 
tively. 

We try to guess that this family corresponds to Vfj. 7/(1^2, W3, W7). Elkies write 
a coordinate x on Vq^/{w2, W3) for which the action of w-j is given by [5, (40)] 
wt(x) = (116 — 9x)/(9 + 20a;). The two fixed points of this action arc —29/10 and 
2. Thus, a coordinate on the quotient is provided by 

The map from Vq^/ (u>2, W3) to Vg* is given by [5, (39)] by 

- (25 + 4a; + 4a; 2 ) (2 - 12 x + 3 a; 2 - 2 a; 3 ) 2 _ , (2a; 2 - x + 8) 4 



(8-4) C = 1/(1 - *) = 



108 (37 - 8 a; + 7 a; 2 ) 108(7a; 2 - 8.t + 37) 

This looks nicer with £ 

108 (37- 8a; + 7x 2 ) 
(8-a; + 2a; 2 ) 4 

The preimage of £ = is 6 times 00 plus the two roots of 7x 2 — 8a; + 37. The 
preimage of £ = 00 is 4 times each of the roots of 2a; 2 — x + 8. The preimage of 
£=1, ort = 0, is 2 times each of the roots of 2 — 12a; + 3a; 2 — 2a; 3 plus each of the 
roots of 4a; 2 + 4a; + 25. Thus the elliptic points of ^6,7/(^2,^3) are the two roots 
of 7a; 2 — 8a; + 37 with index 6 and the roots of 4a; 2 + 4a; + 25 with index 2. two 
pairs of points are interchanged by w?. The elliptic points of order 6 are mapped to 
y = —243/100 while those of order 2 are mapped to y = —24/25. In addition we get 
elliptic points at the ramification points of the covering at y — and y = 00, both 
of order 2. Now, if we guessed correctly, there would be a Mobius transformation 
sending the 4 elliptic points to the 4 singular points of the elliptic surface, sending 
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y = —243/100 to A = . It is easy to check that the unique transformation of this 
type is 

486 + 200y 
~ -216- 225y ' 

Composing with (8.3) we get 

-8 (37 - 8a; + 7a; 2 ) 
A ~~ 9 (25 + 4a; + 4a; 2 ) ' 

This is confirmed by the cr-invariants. Our guess can be proved rigorously [3, 
Subsection 8.5]. 

9. Software 

All the relevant computations for this work can be downloaded from http : 
//www. math. bgu.ac.il/~bessera/picard-fucris/. They are in the form of a 
MATHEMATICA notebook. The notebook is self explanatory. It loads the fol- 
lowing files: 

• pf.m - main file containing all the algorithms 

• data.m - file contains the equations for the elliptic fibrations on Hcrfurtner's 
list 

• elkiesdata.m - contains the differential equations obtained by Elkies. 
The relevant functions contained in the file pf.m 

• picfucs function - computes the Picard-Fuchs equation for an elliptic sur- 
face. This is just a translation into Mathematica of the Maple script 
by Beukers, which may be found at http://www.staff.science.uu.nl/ 
~beukel06/picf uchs .maple 

• Twistpf function - computes the Picard-Fuchs equation for the family of 
twists given the equation for the original elliptic fibration. 

• SigChVar function - makes a change of variable for the sigma invariant. 

• DRes function - computes the residue of a quadratic differential. 
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